Phase diagram of the Hubbard model on honeycomb lattice 
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In this paper we generalized the slave-particle technique to study the phase diagram of the Hub- 
bard model on honeycomb lattice which may contain charge fluctuations. For large U, we have 
antiferromagnetic order phase. As we decrease U below U C 2 — St, the system undergoes a first order 
phase transition into a gapped spin rotation invariant phase. Under a semiclassical approximation 
of the slave-particle approach, we find that such phase breaks the translation symmetry, the parity 
and the lattice rotation symmetry. However, beyond the semiclassical approximation, a Z2 spin 
liquid that does not break any lattice symmetry is also possible. 



Introduction. — Hubbard modelfH, 0] is believed to 
describe the physics of manystrongly correlated sys- 
tems e.g., Mott insulator d, |3] and high temperature 
superconductors 0-0] • It is the simplest model one can 
write capturing the strong correlation physics. So far 
many theoreticalfl, [t| and numerical techniques [Tol - fl^ l 
have been developed to study this model. Among them 
is the slave particle [l3l - [l6j theory which was motivated 
by the RVB state first introduced by P.W. Anderson (l7j. 
One of the interesting phases that have been studied and 
is strongly supported b y t he slave particle approach is 
the Z2 spin liquid phase |18l - [2fj| ] which does not show any 
long range order down to zero temperature. Unfortu- 
nately this phase has not been experimentally verified 
but recently Meng et a/[10] have studied the Hubbard 
model on honeycomb lattice at half filling using the quan- 
tum Monte Carlo (QMC) method and have reported the 
existence of a spin liquid phase for a range of U/t. For- 
tunately QMC does not have sign problem on bipartite 
lattices at half filling so we can trust its results. For 
small U-limit they have reported the semi-metallic phase. 
At U c i ~ 3.5i they have seen a phase transition to the 
spin liquid with nonzero spin excitation gap (gapped spin 
liquid). At U c \ the charge gap opens up and therefore 
this transition point is associated with the Mott metal- 
insulator transition. For a larger value of U C 2 ~ 4.3i 
they have obtained the anti-Ferromagnetic (AF) order in 
which the charge gap is still nonzero but the spin excita- 
tion is the gapless Goldstone mode. 

In this paper, we generalized the slave-particle 
method @, [H| to capture charge fluctuations [HI, [2l[ to 
study the Hubbard model on honeycomb lattice, we have 
obtained a similar phase diagram (see Fig. [1] and [5]) as 
in [To| but with different numerical values for U c i and 
U C 2- We obtained a superconducting phase (instead of 
the semi-metal phase) for small U/t and a AF phase for 
large U/t. At the meanfield level, our phase between 
U c i and U C 2 is a spin liquid with finite charge/spin gap 
that do not break any symmetry. However, the meanfield 
state is unstable. Under a semiclassical approximation, 
we show that phase between U c i and U C 2 is a charge/spin 
gapped state that breaks translation and lattice rotation 
symmetry but not spin rotation symmetry. On the other 
hand, in the presence of the second neighbor hopping, the 
meanfield state may become a Z2 spin liquid state [l8l- 
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FIG. 1. Spin excitation gap (blue line) and charge excitation 
gap (red dots) 



l20( that does not break any symmetry and has finite 
charge/spin gaps. All phase transitions are first order 
which agrees with experiments [4]. 

We would like to point out that the slave-rotor method 
is the other method to include charge fluctuations. [2lj 
which give rise to a nodal spin liquid between 1.68t < 
U < 1.74t. The slave-rotor method is more reliable for 
small U/t and gives rise to the correct semi-metal phase. 
Our method is quite unreliable at small U ft and gives 
rise to a (wrong) superconducting state. 

In the large U/t limit, the Hubbard model can be 
approximated by the Heisenberg model and we expect 
strong AF order in it. This model has been extensively 
studied by different methods (22| - [28{ . Here we use a dif- 
ferent approach to study the antiferromagnetic phase. It 
is shown that the spin/chrage gapped phase has an in- 
stability towards antiferromagnetism. 

The Hubbard model is defined as the following: 

H = Uj2^n iti - t J2 C )« C i° + h - c - (!) 

1 (i,j),<r 

Here means site j is one of the nearest 

neighbors of site i. We know that Hilbert space 
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FIG. 2. Staggered magnetization, m (( — ) l S z (i)), as a func- 
tion of ~. There is a phase transition to the antiferromagnetic 
order at U/t = 3. 



of Hubbard Hamiltonian has four sates per site. 

' lt)i ) \Di > It4-)i- Let's name each state as fol- 
lows: IholonL = h\ |vac). = |0) i , |spin up spinon) ^ = 
fi,t l vac )i = lt) s : i |spin down spinon) . = ft^ |vac) i = 
, Idoublon),; = d\ (vac), = in which |vac) is the 
vacuum, an unphysical state which contains no slave par- 
ticles even holons. Using this picture we can rewrite the 
electron creation operator as: C\ a — j} a hi + a d\f^_ a = 



hi d. 



f 



It should be mentioned that the physical Hilbert space 
contains only four states: empty state (holon), one elec- 
tron (spinon) and two electrons (doublon) on each site. 
So we always have one and only one slave particle on each 
site. So we conclude that we should put the local con- 

if = 1, to get rid of redundant 



straint n? 



states. This is the physical constraint which should be 
satisfied on every site. We could also obtain this result 
by noting that the electron operators are fermion and 
should satisfy the anticommutation relations. From the 
definition of C[ a it is obvious that it is invariant under 
U(l) gauge transformation (We require hi and di to re- 
main bosonic operators i.e., preserve their statistics after 
transformation, otherwise we would have SU(2) gauge in- 
variance. However at U = oo we have only fermions and 
only in that case we have SU(2) gauge symmetry). It is 

I 



worth noting that all the slave particles carry the same 
charge under the internal U(l) gauge. Since the con- 
straint as well as the Hubbard Hamiltonian are gauge 
invariant, so is the action of the Hubbard model. 

In terms of new slave particles, the Hubbard Hamilto- 
nian can be written as: 



t 



(id) 



A/t A b • 



(2) 



notations x 4 3 



f 



A 



r _ 



In which we have used these 
Eit ' Xi.j = hi^j — ^dj . — ( , 

Y,a a f-<y,ifi,° ' = dihj + hidj . To implement the 
constraint we appeal to the path integral and the La- 
grange multiplier methods. 



S = J DpDfDh^DhDd^Dd ■ 



-JdrL 



T - f f — f 

U Ji,<7 Q T Jl, 



d]-^-di 
or 



d 



9* = flhl + flki + h\hi + d\di - 1 



(3) 



The above motivates us to define the effective Hamil- 
tonian as H e ff = H + i^iXigi- Now by using the 
Hubbard-Stratonovic transformation we can decouple 
spinous from [hard-core] bosons at the mean field level. 
To do so we just replace Xi,j an d other operators with 
their average. For translational invariant systems we 
can assume: (Xi,j) = (Xi-j) an d so on - From now 
on x stands for the average of x operators and so on. 
Moreover (iXi) = An. By these assumptions we can ob- 
tain unknown parameters in the effective Hamiltonian 
from self-consistency equations. Now, let us focus on 
the effective Hamiltonian of bosons. As long as A/ is 
nonzero, the pairing between holons and doublons is 
nonzero, and they form bound state. Using the Bo- 
goliubov transformation we can show that the ground- 
state wave-function of bosons is a paired state which 
is completely symmetric between holons and doublons. 
Therefore, as long as this state represents the ground, 



we have (h\ A h k 
Xb = (h\ A hi 



l k.A^ 



dk Bh and as a result: 



d] A d^ B 



l i A n j-B ~ u i A a j,B 1 — 0. Spinons cannot hop 
in this case and the system is insulator. Self-consistent 
equations show that x/ — as well and therefore the 
following Hamiltonians describe the low energy theory of 
this phase: 



H f' B = [ fl,A,t f- k < B 

k 

[ d l,A h -k,i 



H r = 



-A 



-t^k 

+ A 



U - A - tA{ 



-An 



J~k : B,l 

dk,A 
hl k ,B 



(4) 
(5) 



where A / - b ( k\ = J2s &g b e*^ and <5 connects two nearest neighbors. We have similar equations for H 



B.A 
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B A 

and H b ' . Using the Bogoliubov transformation we can 
diagonalize the above Hamiltonians. The energy eigen- 
values for spinons are E[ = \J\l + (tA b k ) 2 . For bosonic 

part we obtain: E±' k = + ± f + \J ( £ ^ a ) 2 - {t^kf ■ 

At half filling, in order to excite a charge we need to an- 
nihilate two spinons and create a pair of holon-doublon. 
So we can define the charge excitation gap as the sum of 
the excitation energy of a holon and a doublon. When 
the charge gap is nonzero then the paired holon-doublon 
state is stable because exciting quasi-particles on top of 
this state costs energy. For this state, the charge gap 

is E« = mmE+' k +minE^' k = 2^/ ' {^^-f - (3tA/) 2 . 
Therefore, as long as U — 2A > 6tA/, charge gap is 
finite and we are in the insulating phase. 

On the other hand, when the charge gap closes, the 
paired holon-doublon state becomes unstable and free 
holons and doublons proliferate. In this state, doublons 
and holons condense independently (single boson conden- 
sation) such that {di^A) — — an d {h^A) — (K,b), 
and therefore Xb — 2 (hi. a) 2 = 2n/i 0, therefore spinons 
can hop freely and the ground state is no longer an in- 
sulator. Since doublons condense at sublattice A and B 
with opposite signs, we show that = and as a result 
Af = 0. we relate the onset of single boson condensation, 
i.e. the critical point below which charge gap closes, to 
the Mott transition. It should be mentioned that Xf,b 
as well as Abj jump at this point, so we obtain a first 
order phase transition in this way, which is consistent 
with experiments. We like to point out that since d)h 
operator that carries 2e electric charge, condenses in this 
state, we indeed obtain a superconducting state instead 
of a semi-metallic phase. 

Phase diagram. — In the following sections we discuss 
the three phases that we have obtained from the slave 
particle method. 

Superconducting phase. — Now let us approach the 
Mott transition point from below i.e. from supercon- 
ducting side. In this phase both Af and A& are zero 
and therefore the charge excitation gap as well as the 
spin excitation gap vanishes. Gapless charge excitation 
implies : muxE^ k +mmE^ k — U — 2A — 6\f = 0. This 
condition can be satisfied up to U c \ = 2A + 6txf = 2.2t. 
At this point the Mott transition happens. In terms of 
physical electrons, we obtain an s-wave superconducting 
state with gapless charge and spin excitations. The pair- 
ing order parameter changes sign under parity and 60 
degrees rotation and transforms trivially under all other 
symmetry transformations. It should be mentioned that 
at small U limit, the Bose gas of holons and doublons 
becomes very dense and there is strong interaction be- 
tween them. So the mean-field results are unreliable in 
this regime and the superconducting state is a fake re- 
sult. However, our method captures two important right 
features of the system below the phase transition, be- 
cause we obtain zero spin excitation energy as well as 
zero charge excitation energy. 



Charge/spin gapped phase. — For U > U& we have 
Xb = 0. So the quasi-particle weight of spinons are 
zero and they cannot hop since for any i and j arbi- 
trary sites: (^fj a fi ja ^ = 0. Therefore this state is like a 
superconductor with infinite carrier's mass m ~ — > 
co. Now let us find U c \. To do so we assume that: 
Af,b (j?J — A/^. So we have Af.b ffcj = Af^r/(k), where 

rj(k) — e lky + 2e~ l ^~ cos^k x and therefore the energy 

spectrum of spinons and bosons are \J \ 2 + jtAj (k)\ 2 and 

±U/2 + \J{U/2 - A) 2 - \tA f (k)\ 2 respectively. From the 
energy dispersion of bosons, one can read that the charge 
gap closes when U c i — 2 A + 6iA/. Our numerical re- 
sults show that near the phase transition, Af ~ .5 and 
A ~ — .4 and the Mott transition occurs at U c i/t = 2.2. 
For large U/t limit: A f -> .53 , A b - £ , A - (^) 3 Ln^ 
and rib ~ A^ ~ {jj) 2 ■ It is clear from the energy spec- 
trum of spinons that in the spin liquid phase, there is a 
gap in their spectrum equal to: E^ — |A|. Note that in 
the spin-charge separation picture, the physics of spin is 
determined by that of spinons. Therefore the spin exci- 
tation gap is also Eg = |A|. 

Now let us focus on the gauge theory of this phase. In 
this phase the effective action of spinons is of the follow- 
ing forms: 

= ^ X/ fi,T,ofi,T,a 

-t AM'-./W: ,.,,/:/;. „ • h,: (6) 

<i, - ?>,(T,r 

Now if we transform operators as: fi } A t a & la fi,A,a 
and fi^B.a s la fi,B,a for any arbitrary phase a, i.e. 
assuming a staggered global gauge transformation, then 
the effective Hamiltonian does not change. Therefore the 
invariant gauge group (IGG) of the Hamiltonian is the 
staggered U(l). The reason is that there is no hopping 
term due to the nonzero charge gap and the gauge trans- 
formation of two neighboring sites have opposite phases, 
the total phase change of the pairing term becomes zero 
and therefore gauge fluctuations are described by stag- 
gered compact U(l) instead of compact U(l) gauge the- 
ory. This is equivalent to assuming positive unit charge 
on sublattice A and negative unit charge on sublattice B 
for slave particles under the internal gauge transforma- 
tion. 

So, at mean field level, the charge/spin gapped phase 
has a neutral spinless U(l) gapless mode as its only low 
energy excitations. However, it is well known that 17(1) 
theory in 2+1D is confined due to instanton effects. So 
let us assume that the U(l) fluctuations are weak and 
use the semiclassical approach to study the U(l) confined 
phase where the U(l) mode is gapped. We find[29| that 
these instanton operators, e 10 (in the dual XY model), 
carry a non-trivial crystal momentum. Also, under 60 de- 
gree lattice rotation and parity, an instanton is changed 



4 



to an anti-instanton, e 10 — > e~ %6 '. The instantons carry 
trivial quantum numbers for other symmetries. However, 
a triple instanton operator cos(3#) carries trivial quan- 
tum numbers for all symmetries. This allows us to con- 
clude that the neutral spinless U(l) mode is described 
by L — j^(dO) 2 + If cos (3$). In the semiclassical limit 

(the small g limit), (e* 6 ') ^ and we obtain a phase that 
breaks the translation, the parity and the 60° rotation 
symmetries, but not spin rotation symmetry. 

We like to point out that in the presence of second 
neighbor hopping in the Hubbard model the charge/spin 
gapped phase can be spin liquid that do not break trans- 
lation, parity, 60 degree lattice rotation, and spin rota- 
tion symmetries. It is because we can break the stag- 
gered compact U(l) gauge symmetry down to a Z 2 one 
by Anderson-Higgs mechanism. If we add second neigh- 
bor hopping to the Hubbard model, within slave par- 
ticle approach, this term generates pairing terms of the 
form fJ T a fj T _ CT , i.e. it induces the same sublattice pair- 
ing and the Hamiltonian is no longer invariant under the 
staggered global U(l) gauge transformation. In this case 
the staggered compact U(l) gauge symmetry is broken 
down to a Z 2 gauge symmetry. The U(l) gauge fluctua- 
tions are gapped and thus our mean filed state is stable 
and we can trust our meanfield results. Therefore we 
obtain a spin liquid phase. 

Antiferromagnetic phase. — In this part we show that 
the charge/spin gapped phase is unstable towards an- 
tiferromagnetic order above U C 2 — 3t. To obtain Neel 

order in the t-J model we simply assume that (S z a 



^S z .bJ — m. But how can one implement this idea 
in the Hubbard model within slave particle approach? 
In the Neel order phase, translation symmetry is bro- 
ken and there is an asymmetric situation between sub- 
lattice A and B. For example we can obtain a antifer- 
romagnetic phase by assuming Axj = (fj,B,ifi,A,t) 
{fj,A,lfi,B,t) = ^2./- This assumption simply means 
that the chance of finding a spin-up spinon on sublat- 
tice A and another spin-down spinon on sublattice B 
is more than finding the opposite one, so this method 
introduces staggered sublattice magnetization and leads 
to the Neel order. If there is a Neel order in the 
system then the chance of creating one holon-doublon 
pair from annihilating a spin-up spinon on sublattice A 
and a spin-down spinon on sublattice B is more than 
the other process. Therefore the excitation energy of 
spinons for up-spin on A and down-spin on B is E^{k) = 



y/X 2 + |iAi.fc(fc)| 2 while for down-spin on A and up-spin 
on B is Ej(k) = y/X 2 + |tA 2 , b (fc)| 2 . On the other hand, 
since we are not interested in CDW, we need a sym- 
metric situation between sublattices A and B for the 
charge sector. So the energy excitation of bosons is 

E b (k) = \J{U/2 - X) 2 - |tA/(fc)| 2 . Using these assump- 
tions we lead to the following self-consistency equations: 



Ai,/ 



N. 



t x ^ |77(fc)| 2 Ai, 6 



t v |77(fc)| 2 A 2 , b 
2J N s ^ E 2J (k) 

A f = A 1J + A 2J 
A hb A hf + A 2 , b A 2J 



t v \y(k)\ 2 A f 
N s 4" E b (k) 



(7) 

(8) 

(9) 
(10) 



By solving the above equations we find that above 
U C 2 = 3t, m 7^ 0. So we conclude that for U > U c2 
we obtain AF order. It is interesting that in this phase, 
the gap of spinons is very small and negligible (for ex- 
ample at U=4, it is —2 x 10~ 7 ). So in this phase we can 
assume that spinons are massless quasi-particles. 

In conclusion, we have used a generalized slave parti- 
cle method to derive the phase diagram of the Hubbard 
model at half filling on the honeycomb lattice. Within 
the mean field approximation we can decouple fcrmions 
from bosons to achieve the effective Hamiltonian that de- 
scribes the low energy physics of the system. The physics 
of the Mott transition is discussed and it turns out to be 
a first order phase transition. It is shown that the phase 
transition occurs when the charge gap opens up. Above 
the critical point, within meanfield theory we obtain a 
spin liquid phase. But after including gauge fluctuations 
of the emergent spin liquid and investigating the instan- 
ton effect, we argue that this phase is unstable and we 
finally obtain a spin/charge gapped phase that breaks 
the translation symmetry. For large U limit, a new ap- 
proach to study antiferromagnetic phase within the slave 
particle picture has been developed. It is shown that 
the gapped spin liquid phase has an instability towards 
antiferromagnetism. 
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